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Applying multicanonical simulations we investigated folding properties of off-lattice het- 
eropolymers employing a mesoscopic hydrophobic-polar model. We study for vaiious se- 
quences folding channels in the free-energy landscape by comparing the equihbrium confor- 
mations with the folded state in terms of an angular overlap parameter Although all investi- 
gated heteropolymer sequences contain the same content of hydrophobic and polar monomers, 
our analysis of the folding channels reveals a variety of characteristic folding behaviors known 
from realistic peptides. 



1 Introduction 

The identification of folding channels is one of the key tasks of protein folding studies. 
While secondary structures depend on atomistic details (such as, in particular, hydrogen 
bonding), tertiary structure formation should exhibit a certain degree of universality. This 
suggests that coarse-grained models might capture the main characteristics on mesoscopic 
scales of this stage of the folding process. In this note we report a computer simulation 
study that tests this idea by employing the off-lattice hydrophobic-polar AB model. 



2 Model and Method 

In the AB model' a heteropolymer or coarse-grained peptide is described as a chain of 
hydrophobic (A) and hydrophilic (B) monomers whose energy is obtained from specific 
Lennard- Jones potentials between all nonbonded pairs. Thereby an aqueous environment 
is modeled implicitly by the energetically favored A-A contacts that lead to the formation 
of a hydrophobic core at low temperatures. A smaller energy contribution arises from the 
local bending of the chain. 

In the Monte Carlo simulations we applied the multicanonical technique^ where an 
additional weight function leads to a flat distribution in energy space. This allows the 
reweighting of the data to any desired temperature with equally high accuracy. In order 
to identify folding channels we used the angular overlap parameter Q introduced in Ref. 3 
to compare equilibrium conformations at temperature T with the previously determined 
ground-state structure. Two conformations are equal if Q = 1. The temperature-dependent 
probability distributions of the total energy E and Q, Pt{E, Q), then allow the analysis of 
the chain's folding behavior (for more details, see Ref. 4). 



3 Results 



Folding through intermediates: For the sequence A4BA2BABA2B2A3BA2 the conforma- 
tions at higher temperatures do not exhibit significant similarities {Q k, 0.7) with the 
ground-state conformation depicted in Fig. 1. As the distributions Pt{E,Q) in Fig. 2 
show, the main branch sHghtly moves to higher Q-values with decreasing temperature un- 
til it splits below T « 0.2. Near T w 0.1, the population of intermediate conformations 
has increased and coexists with denatured states. Approaching T « 0.05, the intermediate 
states dominate the canonical ensemble. The probability for denatured conformations is 
reduced, but the onset of an occupation of states with similarities to the global-energy min- 
imum (GEM) is clearly visible. At T = 0.02, the majority of conformations has a large 
overlap with the ground state and the heteropolymer has folded into its native state. 




Figure 1 . Ground-state conformation of A4BA2BABA2B2A3BA2. 



T = 0.2 T = 0.1 




Figure 2. Canonical distributions Pt{E,Q) of energy E and overlap Q with the GEM for the sequence 
A4BA2BABA2B2A3BA2. 
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Two-State folding and metastability: A typical two-state folding scenario can be observed 
for the sequence BAgBA4BA2BA2B2, where again Q « 0.7 at high temperatures indicates 
that the heteropolymer is in a random state that possesses no similarities with the GEM 
conformation. Around T = 0.1 the canonical ensemble divides and a state with Q > 0.9 
is occupied in addition to the disordered conformations with Q < 0.8. With decreasing 
temperature the newly formed branch dominates the ensemble and approaches Q = 1.0 
when T ^ 0. 

An example for metastabiUty is provided by a third heteropolymer with the sequence 
A4B2A4BA2BA3B2, which exhibits in the low-temperature phase a glassy behavior. While 
at high temperature similar to the other cases all conformations are distributed around 
Q » 0.7, the main folding channel does not lead to a single ground-state conformation. 
Instead two rival conformations (Q » 1 and Q i=s 0.75) can be found also at temperatures 
below T = 0.01. 

4 Summary 

Numerical investigations of medium size or large proteins by means of Monte Carlo sim- 
ulations are very difficult since complex energy functions demand high computational ef- 
forts. On the other hand, problems of comparable complexity with equal or higher number 
of degrees of freedom can be handled using simpUfied coarse-grained models Uke the AB 
model. Although the model is relatively simple and the investigated sequences are only 
permutations of each other (and were not designed especially), three different kinds of 
folding could be observed. Since all of them - folding through intermediates, two-state 
folding and metastability - are also known from real peptides, the AB model seems to 
resemble general characteristics of protein folding. Therefore we believe that further re- 
search on this model offers the possibihty to gain quahtative insights in tertiary folding 
where microscopic details are expected to be of less importance. 
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